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FOREWORD 


A method to determine the modes of dynamical systems is 
described. The inputs and outputs of a system are Fourier 
transformed and av'eraged to reduce the error level. An 
instrumental variable method that estimates modal parameters from 
multiple correlations between responses of single input, multiple 
output systems is applied to estimate aircraft, spacecraft, and 
off-shore platform modal parameters. 
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I. INTRODUCTION 


A major task in the development of structures (ground 
structures, aircraft and spacecraft) and their control systems is 
the identification and demonstration of satisfactory dynamic 
characteristics . Each specific structural system, however, has 
unique characteristics which require special analytical 
treatment. For the development of aircraft, it is not only 
necessary to determine its flutter boundary or establish that the 
aircraft is flutter-free throughout its operating .nvelope; it is 
also desirable to obtain detailed information about dynamic 
characteristics such as aircraft gust response and provide data 
for the verification of mathematical models of the aircraft. The 
determination of the dynamic characteristics of off-shore 
platforms is essential for the verification of their safety in 
different sea states. The modes of large space structures are 
very close to each other, and it is necessary to identify them 
very accurately not to excite them inadvertently by the control 
system [1,2). The most important parameters are generally those 
defining the frequency and damping of the low-damped modes of the 
system. A number of surveys show the development of stability 
testing of aircraft [3-7]. These procedures typically involve 
fast sine-sweep excitations and analysis of the response. In 
Refs. 8 and 9, methods for the identification of frequencies and 
mode shapes of off-shore platforms have been developed. These 
methods are based on the maximum entropy spectral analysis method 
[ 10 ]. 

Testing of the dynamic characteristics of systems involve 
three tasks: design of an excitation input that excites all 
significant modes, measuring the response of the system, and then 
extracting the modal parameters from the response. In Ref. 11, it 
was shown that the best accuracy in modal parameters can be 
obtained for a given amount of data by exciting the system by an 


external input in the required frequency band and measuring the 
system transfer function. 

Testing methods, such as fast sine-sweep measurement of the 
transfer funtion and measuring the transient decay after an 
impulse input, rely on a high signal-to-noise ratio and obtain an 
accurate response measurement in a laboratory environment [11]. 

The errors of these measurements have only small correlations and 
can be considered to be white. In-flight tests cf the measured 
accelerations and moments of aircraft are not only the response to 
excitation inputs of the control surfaces, but also to gust 
inputs. The gust responses are necessarily non-Gaussian 
distributed and highly correlated. In Refs. 15 and 14, it was 
shown that the least squares estimates are very sensitive, to the 
Gaussian assumption and the mean estimate is much more robust. 
Reference 11 shows an example, where correlating the response and 
a measured external .input, as well as averaging reduces the test 
time, for a given level of accuracy by a factor of 10 to 15. 
Reference 15 combines the robustness of the single-input/ 
single-output method of Ref. 11 with the multi-input/multi-output 
method of Ref. 12 to design a robust and highly efficient method 
for flutter testing. The robustness of this algorithm is due to 
averaging of the autospectra of different measurements. 

References 11 and 16 show that the autospectra, as well as the 
cross spectra of measured outputs contain information about the 
modal parameters of the test structure. 

This report is concerned with estimating modal parameters by 
combining autospectra and cross-spectra of responses for different 
time delays. In the following sections, a frequency domain 
derivation of an instrumental variable algorithm is shown. Four 
examples of its application to modal parameter estimates of 
aircraft, an off-shore platform, and a simulated space structure 
are presented. 


II. INSTRUMENTAL VARIABLES ALGORITHM 
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The input/output relation of a system can be described to 
first order by 

y » X0 + V (2.1) 

where y and X are measured data, 6 are the system parameters 
to be identified, and v is some zero mean noise term which is 
assumed to be white. The least squares estimate of e is 

®LS “ (x'^X)'^x'^y (2.2) 


Ele^^gJ » e + (X’^X)'^ E(Xvl. 

The least squares estimate, therefore, has a bias of the 
measurments X and the noise v are correlated. In the 
instrumental variables method Eq. (2.1) is premultiplied by a 
matrix [Zjrank 2 = rank X] that is not correlated with v. 

Hence, the instrumental variables estimate of e , 

» (2^X)“h'^y, 

E{9jy] » e (2.3) 

is unbiased. In the present case, we use for instrumental 
variables measurements of a past time. Choosing a time delay, 

T , sufficiently large these measurements are not correlated with 
the noise of the current ones. 

The input/output relationship of a linear t ime- invariant 
dynamic system can be modeled by a set of linear constant 
coefficient differential equations (Appendix A). In the frequency 
domain these equations may be written as 


. y„ Kjo,) 
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are the coefficients of the characteristic polynomials for poles 
and zeros, respectively 112), and and are the 

Fourier transforms of output and input, [y(t), u(t)l 0 < t< T). 

In Refs. 12 and 15, Eq. (2.4) was premultiplied by 
IC-joj) ... IJ Yq e-* 


to get a set of independent equations. Due to the o)>dependent 
phase of Y^*eJ‘"^, this causes some frequencies to be 
weighted more heavily than others and is the reason for some of 
the difficulties in estimating ARMA parameters for finite data 
sets. By premultiplying Eq. (2.4) instead by the instrumental 
variable, 

Y-t " F[y(t-T)|- X 1 t £ T-t], 


the undesirable frequency weighting is avoided. This follows from 
the fact that 

ElYoI - E(Y.^ ej“^] 

for a stationary system. Premultiplying Eq. (2.4) with 

gives an additional independent equation for each of the p 
outputs. The resulting equations are shown in Appendix A. 

There are several methods available to solve the instrumental 
variables ARMA equations; some of which include least squares, 
averaging, maximum likelihood and prediction error methods. The 
maximum likelihood method and the prediction error method 
generally give the best fit but their complexity and computation 
requirements make them unsuited for on-line applications. The 
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least squares and averaging methods are fast and well suited for 
on-line applications IlS). 


The Iv'ast squares method of Appendix B estimates ARMA models 
that fit the measured spectra very well, but is very sensitive to 
non-Gaussian noise and spurious peaks in the spectra. The 
averaging method of Appendix C reduces the parameter uncertainty 
for small data sets, but increases the fit error. A flowchart of 
an on-line application of this algorithm is shown in Figure 2.1. 


INPUT & CHfTPUT 
TIME HISTORIES 



Figure 2.1 Schematic Flowchart of the On-Line 
Dynamic Modal Estimation Algorithm 
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III. MODEL ORDER DETERMINATION 


In real-time estimation of dynamic stability, the number of 
important modes is generally not known. Therefore, it is often 
necessary to estimate the number of significant modes from the 
input/output measurements. Akaike [16] determines the model order 
of a system by finding the minimal state-space realization from 
the associated Hankel matrix. Tse and Weinert [17] find the 
minimum order by successively increasing the order of a state- 
space realization until the determinant of the associated 
prediction of the output covariance shows a sharp drop kn 
magnitude. Rissanen [22] used an entropy argument in a 
theoretical discussion on model structure identification. Young 
and Jakeman [23] use the largest eigenvalue and the weighted trace 
of the parameter error covariance to determine model orders for 
systems with deterministic inputs. These methods show good 
results if the disturbances of a system are white. Due to the 
non-white nature of gust disturbances in aircraft testing it is 
extremely difficult to estimate the modal order by the above 
methods . 

A fundamentally different approach to the problem of 
estimating the order of the model based on input/output 
measurements is proposed here. The basic idea is as follows. If 
the specified model order is lower than the true model order, the 
identification algorithm in this chapter will give incorrect 
estimates of the poles and zeros. If the specified model order is 
the same as the true model order, correct poles and zeros will be 
identified. However, if the specified model order is higher than 
the true model order, there will be many superfluous poles and 
zeros in addition to the true poles and zeros. The superfluous 
poles must approximately cancel out the superfluous zeros. 
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Therefore, if we specify a maximal model 
this chapter, the true model order could 
transfer function is simplified first by 
zeros which are very close to each other 
and zeros outside the frequency range of 


order in the algorithm of 
be obtained after the 
cancelling the poles and 
and, secondly, the poles 
interest . 
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IV. EXAMPLES OF INSTRUMENTAL VARIABLE MODAL 
ESTIMATION APPLICATIONS 


Instrumental variable modal estimation can be applied to many 
different model structure and parameter estimation problems. The 
examples shown in this chapter were chosen for their diversity and 
availability of data. The flight test data of the B-S2, supplied 
by the Air Force Flight Dynamics Laboratory (AFFDL), are of high 
quality and are well-suited for dynamic structure identification. 
The flight test data of the Boeing 747 SCA - Space Shuttle 
combination include missed data points. The off-shore platform 
data include no measurements of the inputs and therefore require a 
significant amount of averaging to reduce biases in the 
estimates. The space structure data are simulation data with 
white noise and without any nonlinearities. 

4.1 APPLICATION OF THE MODAL ESTIMATION ALGORITHM TO B-52 FLIGHT 

TEST DATA 

Interaction of structural flexibiity with flight modes has 
been shown to be very important for large and flexible aircraft. 
The low drag and highly maneuverable design of modern aircraft 
make them even more flexible and put stringent requirements on the 
control system design. The control system not only has to provide 
aircraft stability, but also to control elastic modes to relieve 
stresses and pre snt flutter. These control configured vehicles 
(CCV) concepts of modern aircraft can only be realized if the 
aircraft dynamics in general and particularly the aeroelastic 
behavior is known sufficiently accurately. The AFFDL conducted 
flight tests with a B-52 CCV specifically designed to provide data 
for dynamic model and parameter estimation (18]. 

Accelerometers and gyros were located at different parts of 
the aircraft (Figure 4.1 and Table 4.1) to measure accelerations 
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Figure 4.1 Sensor Locations of B*52 


CCV 



Figure 4.2 Locations of Bending Moment Meaurements 

on B-52 CCV 








Table 4.1 

Recorded Flight Test Measurements 
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ITEM 

iXA' ! 

1 

V«rt1C4l EtiMIng Montflt 

SS761! 

2 

Vtrticti NnOing Nontnt 

BS1222 

i 

^trtiC4l tinging McNivnt 

BS1412 

4 

L4tirt1 Binding MoiMnt 

BS760 

S 

Littrti Binding Moniint 

851222 

4 

Latirti Binding HoMtnt 

BS1412 

7 

Vartlcal Binding Mafwnt 

IWS222 

8 

Vartlcal Banding Mominc 

IW5820 

9 

Vtrtlcal Binding Monint 

IWS974 

10 

ChordHl*i Binding Motnint 

lilS222 

n 

Cnordwiti Binding Menwnt 

l'4SB20 

12 

Chor*dMlsi Banding Nenanc 

IW5974 

13 

Vartlcal Banding Nomant 

AUS222 

U 

Vartlcal Sanding Monant 

8WS820 

15 

Vartlcal Banding Momant 

8W5974 

18 

ChordMiti Banding Momant 

BU5222 

17 

:*^rdw1sa Binding Ho*»nt 

AWS820 

18 

''•‘dwlia Banding **0inint 

RWS974 

19 

leal Banding Moatant 

IHT8156 

20 

Vartlcal Banding Monant 

AHT81S6 

21 

Vattlcal Sanding Moaant 

PSISS 

22 

aittn itata 

Bsaeo 

23 

toll lUca 

85880 

24 

Yaw Bata 

SS880 

25 

vartlcal Accalaratlon 

85172 

26 

Vartlcal Accalaratlon 

BS860 

27 

Vartlcal Accalaratlon 

351655 

23 

(.ataral Accalaratlon 

35172 

29 

Lataral Accalaratlon 

B5860 

30 

Lataral Accalaratlon 

B51655 

31 

Vartlcal Accalaratlon 

IWSL565 

32 

Vartlcal Accalaraticn 

LMBL925 

33 

vartlcal Accalaratlon 

IWSI 359 

34 

Vartlcal Accalaratlon 

left Ext. Tank, No$a 

35 

Vartlcal Accalaratlon 

AUH56S 

36 

Vertical Acceleration 

AUB1925 

37 

Vertical Acceleration 

toSI 359 

33 

Vertical Acceleration 

Aignt Ext. Tank, tola 

39 

vertical Acceleration 

left Inboard tocalle 

40 

lataral Accalaratlon 

left Inboard tocalle 

41 

Vertical Acceleration 

left Outboard riacella 

42 

lateral Accalaratlon 

left Outboard tocalle 

43 

Angular Accalaratlon, Y>Ax1s 

LW5540 

44 

Angular Acceleration, Y-Axlt 

US1359 

45 

Angular Acceleration, Y-A*li 

AU51359 











and angular rates for mode identification; strain gages 
attached to measure structural deformations (Figure 4.2 
4.1) to assist in the identification of mode shapes. 






were 

and Table 


For the flight test investigated, the aircraft was flying at 
an altitude of 21,000 feet at 305 KCAS. The excitation was a 
verticil cannard sine swe®’^ 'om 0.5 to 5.0 Hz. Figure 4.3 shows 
response data of this flight, test. In the application ofthe 
instrumental variables modal estimation algorithm (IVME) to these 
data, :he frequency range of interest i'* limited to 
0.1 < oj < 5 Hz and the number of poles and zeros in each transfer 
function are specified to be eight and seven, respectively. The 
minimum time delay for the instrumental variable is. 


T = 2nAtjn = max(poles, zeros). At ■ time between 

measurements 


Transfer function measurements with t, 2t, and 3t where combined 
to improve the estimation accuracy. Discrete frequency noise was 
eliminated from the data set. Figures 4.4 and 4.5 show the 
measured transfer functions between te vertical cannard and 
wing-tip and wing-root, respectively. Elastic modes, as well as 
aeroelastic modes, are clearly excited in the frequency range of 
the si-'“ sweep, as observed in the magnitude plots, as well as the 
phase plots. Table 4.2 compares the theoretical structural modes 
[18] with the identified frequencies and damping ratios. Some of 
the identified modes (4, D, 6, G, 9 H) mask the theoretical 
structural, others (3, C, 7, 8, F) could be low-damped aeroelastic 
modes. 


4.2 APPLICATION OF THE IVME ALGORITHM TO IDENTIFY THE DYNAMICS 
OF THE BOEING 747 - SPACE SHUTTLE COMBINATION 

Boeing has modified the model 747-123 for NASA to carry the 
Space Shuttle Orbiter aircraft "piggy-back" fashion (Figure 4.6). 
The dynamics and aerodynamics of this aircraft combination are 
significantly different from that of the regular B747 and required 
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Table 4 . 2 

Comparison of Theoretical and Identified Modes 
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a new control system design incorporating all the significant 
modes of the B747-Space Shuttle (B747 SCA). Flight tests were 
performed by Dryden Flight Research Center to determine the flight 
envelope of the B747 SCA and to fine tune the control and 
stability system. 

This investigation applies an instrumental variables modal 
estimation algorithm to flight test data of the B747 SCA to 
determine its structural and aeroelastic characteristics. 

For the flight test investigated, the aircraft was flying at 
an altitude of 22,000 feet at M ■ .56. The excitations were 
rudder and elevator sine sweeps. Figures 4.7 and 4.8 show the 
power spectra of the rudder and elevator excitations, 
respectively. These excitations are band-limited and are not 
constant over the range of excitation frequencies. Table 4.3 
shows some of the recorded flight test measurements. For this 
investigation wing-tip acceleration is used to identify the B747 
modes, and the orbiter attachment accelerations to identify the 
modes of the B747-Spice Shuttle combination. Figures 4.9 and 4.10 
show the measured transfer functions for the rudder sweep and the 
elevator sweep, respectively. The measurements show several 
extremas for missed frequencies during the frequency sweep. 

Lacking are good measurements for low frequencies, but some modes 
can be clearly observed. Tables 4.3 and 4.4 show the identified 
antisymmetric and symmetric modes, respectively. Even so, the 
measured transfer functions are relatively noisy; the identified 
modes change only little for different assumptions for poles and 
zeros of the transfer function. These consistent estimates 
demonstrate the robustness of instrumental variables estimates, 
specifically when using multiple correlations for instrumental 
variables of different time delays. 


18 


Table 4.3 

Recorded Flight Test Measurements of B747 SCA 


r 


1 . 

1 

j 

I 


NO. 

location 

UJ 

1 


Pitch Angl* 

2 


Side Slip Angle 

3 


Side Slip Angle 

4 

AR3300A 

Pitch Rate CG Flight Test Gyro 

5 

AR3301A 

Roll Rate CG Flight Test Gyro 

6 

AR3302A 

Yaw Rate CG Fllgni Test Gyro 

7 

AR09110 

Pitch Angle CG Ins No. 1 

3 

AH8S370 

Roll Angle CG Ins. No. 1 

9 

AH95330 

Yaw Angle CG Ins Heading, True 

10 

VTA 

True Air Speed 

11 

QBAR 

Sarometrlc Pressure 

12 

? 

Altitude 

13 

IH1176A 

L Inboard Aileron Pos 

14 

IH1180A 

R Inboard Aileron Pos 

IS 

IH1189A 

L Outboard Elev Pos Act 

16 

:h1192A 

L Inboard Elev Pos Act 

17 

IH1198A 

R Inboard Elev Pos Act 

18 

:h1199A 

R Outboard Elev Pos Act 

19 

IH1207A 

Lower Rudder Pos Act 

20 

IH1211A 

Uooer Rudder Pos 

21 

IH1219A 

No. 12 Spoiler Pos 

22 

IH6*08A 

No. 1 Spoiler Pos 

23 

AA3360A 

Lat CG Accel 

24 

AA'9S13A 

Normal CG Accel 

25 

A;,3632A 

Longitudinal CG Accel 

26 

AA9S14A 

Normal Cockpit Accel 


'1A8332A 

Lat Cockpit Accel 

mSM 

RA1926A 

Lat L Aft Orb Attach ^tng Acce 

29 

RA1927A 

Longitudinal Fwd Orb Attch Ftn 

30 

RA1956A 

Normal Fwd Orb 'Jo Atch Ftng Ac 

31 

RA1957A 

Lat Fwd Orb Attach Ftng Accel 

32 

RA1953A 

Normal L Aft Oro Up Aten Ftng 

33 

RA1961A 

Normal R Aft Orb Jp Atch Ftng 

34 

RA2567A 

Lat L Tip Fin Lower cs Accel 

35 

RA2563A 

Lat L Tip Fin Upper FS Accel 

36 

RA2593A 

Lat R Tip Fin Upper FS Accel 

37 

RA4S06A 

Normal Aft Fus 3S 2300 Accel 

38 

RA6149A 

Lat Tip Vert Tall FS Accel 

39 

RA6151A 

Lat Tip Vert Tall RS Accel 

40 

RA6224A 

Lat Eng 1 Forward Accel 

41 

RA6225A 

Lat Eno 2 Forward Accel 

42 

RA6324A 

Normal Eng 1 Fwd Accel 

43 

RA6325A 

Normal Eng 2 Fwd Accel 

44 

RA6375A 

Normal R Horiz Tail FS Accel 

45 

RA6377A 

Normal P Horiz Tail RS Accel 


RA6379A 

Normal L Horiz Tall *^S Accel 


RA6381A 

Normal L Horiz Tall RS Accel 

43 

RG6384A 

Lat Stab Jack Screw Accel 

49 

^G6386A 

Long L Stab T'p Accel 


RA6392A 

Normal R king Tip RS Accel 


RA6339A 

Norra’ L wing Tip RS Accel 


RA6393A 

Normal L wing Tip FS Accel 


RA9442A 

Lat Aft Fus 3S 2300 Accel 

54 

XX2569 

LH Tip Fin L'NR Long Acc 

55 

XX2570 

LH Tip Fin Upr Long Acc 

56 

XX3361 

Normal Accel 

£7 

IH 53A 

L Outboard Aileron Pos at Ac: 

53 

IHllSlA 

R Outboard Aileron ?os 

59 

Ih671 1 a 

No. 3 Spoiler Pos 

60 

IH6705A 

No. 5 Scoller Pos 

61 

IH7293A 

No. 3 Scoiler Pps 

62 

Inl218A 

No. 13 Scoiler Pos 













I'igurc 4.7 Power Spectrum of Lower Rudder Position of B747 SCA 
(Flight 01-25-77, 12:28:47 - 12:29:22, At =0.1 sec) 






:trum of Right In-Board IMevator Position of B747 
L-25-77, 12:31:30 - 12:33:11, At= 0.1 sec) 







Trr 


T¥7i 5T:T^ 




TiTi — j?:r 

W IRIkD/SCC) 


Trr 


ITT 


^WTT 


ia«. 


1 1 « . 


Figure 4.9a Measured Transfer Function of B747 SCA Be- 
tween Lower Rudder and Lateral Forward 
Orbiter Attachment Fitting Acceleration 
(R-Sweep 01-25-"7, 1024 Points/Record J 
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Measured Transfer Function of B747 SCA Be- 
tween Lower Rudder and Differences of Ncrmal 
Accelerations of Left and Right Aft Upper 
Orbiter Attachment Fittings (R-Sweep 01-2S-7 
1024 Points/Record, Three Averages! 
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Figure 


.9c Measured Transfer Function of B^i" SCA Be- 
tween Lower Rudder and Differences of Wing 
Tip Accelerations (R-Sweep 01-25-7", 1024 
Points/Record, Three Averages) 
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Figure 4.10b Measured Transfer Function of B747 SCA Be- 
tween Elevator and c.g. Pitch Rate (E-Sweep 
01-25-77, 1024 Points/Record, 12 Averages) 
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Table 4.4 


Identified Antisymmetric Modes of the B747 SCA 
l.R-Sweep 01- 25- 77 , 1024 Points/Record) 



Table 4.5 

Identified Symmetric Modes of the B"47 SCA 
(E-Sweep 01-25-77, 1024 Points/Record) 


























4.3 APPLICATION OF THE IVME ALGORITHM TO OFF-SHORE PLATFORM DATA 


The accident of the oil platform off the coast of Norway this 
year shows the importance of having a good understanding of the 
dynamics of off-sho^e platforms and its interaction with the ocean 
environment. It is of particular interest to know the structural 
modes and damping characteristics of the platforms. Measuring and 
identify ng these eigenmodes will assist in determining the safety 
of new platforms at different sea states and points towards the 
need for additional passive or active damping of existing 
platforms . 

In Ref. 8, J.K. VanDiver and R.B. Campbell used maximum 
entropic spectral analysis to identify structural modes of 
off-shore platforms from accelerometer measurements as shown in 
Figure 4.11. The instrumental variables modal estimation 
algorithm is applied here to the same data used in Ref. 8. The 
input to the structure is assumed to be random noise. The 
zero-input variation (A-7) of the algorithm therefore is used for 
this problem in mode identification. Since no phase information 
is available from these data, the accuracy of the mode estimates 
is much reduced and many data sets have to be averaged to achieve 
sufficient certainty in the estimated eigen- f requenc ies and 
damping factors of the structures. 

The resolution of data in the frequency domain increases with 
the observation time. The maximum frequency is proportional to 
the sample rate. The accuracy of mode estimates at low 
frequencies therefore can be improved by increasing the number of 
independent data sets, to be averaged, by jumping over data 
points. If the sample rate is r = 1/at, but the maximum 
frequency of interest is 

n^ = integer (-j^) 

data sets can be averaged. For high-order systems it is necessary 
to run the data sets though a band limited filter to avoid 
alienation due to higher frequencies. 
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Figures 4.12 and 4.13 show the spectra for 2048 data/record, 
28 records averaged and 1024 data/record, 44 records averaged, 
respectively. The spectra in Figure 4.12 have a higher resolution 
than the spectra of Figure 4.13. They are, on the other hand, 
more noisy and give a larger variance in the damping-ratio 
estimates. Table 4.6 compares the estimated eigen-frequencies and 
damping ratios of the off-shore platform with the results of 
VanDi ver/Campbell [8] and the theoretical values [8]. The 
estimated modes determined by IVME are similar to the estimates 
determined by the MEM [8,10]. Both of these estimates are 
significantly different from the theoretical modes of the elastic 
model structure of the platform. This clearly shows the necessity 
of identifying eigen-frequencies and damping ratios of installed 
platforms to determine the safety of the platform at different sea 
states as well as the requirement for passive and/or active 
damping elements. 

4.4 APPLICATION OF THE IVME ALGORITHM TO SPACE-STRUCTURE MODAL 

ESTIMATION 

In the past, the flexibility of spacecraft had a minor impact 
on the mission performance. Analysis was largely confined to load 
analysis during launch and the structural interaction with the 
attitude controller. The spacecraft could, in principle, be 
treated as rigid with flexibility in the appendages. Techniques 
for stabilizing the attitude usually relied on frequency 
separation, where the lowest mode was kept separated from the 
control bandwidth. Large space structures (LSS) will be erected 
in space, free of the constraints on configuration imposed by the 
launch environment. Spacecraft no longer retain the central rigid 
body, but require distributed structural flexibility. LSS have 
more similarities with a large leaf of thin paper than the 
classical point-mass model of previous spacecraft. Closed-loop 
controllers, which apply forces to the structure based on sensed 
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Figure 4.12b Spectrum of Acceleration in y Direction 
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Figure 4.12c Spectrum of Rotational Acceleration 
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Figure 4.13a Spectrum of Acceleration in x Direction; Number of Data Points 
per Record = 1024; Number of Averages - 114; Number of Shifted 
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Figure 4.13c Spectrum of Rotational Acceleration 
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movement, can alter the vibration modes and shapes significantly. 
Active control for LSS is, therefore, not only used for damping 
and attitude control, but also stiffening of the structure and 
shape control. 

The problem of controlling a LSS is difficult because the 
dynamic model is infinite, whereas the physically realizable 
controller must be finite, both in space and the number of 
observations required. The finite dimensional model is generally 
determined by finite element modeling techniques and reduced for 
the control realization. Finite dimensional models generally only 
include the lowest order modes [20] and the modes close to 
excitations (machine vibrations, solar, magnetic, residual 
aerodynamic or gravity forces). As shown in Ref. 19, even for 
this radically reduced order system the poles and zeros are very 
closely placed. Even small misalignments and/or model uncertainty 
of the robust colocated sensor-actuator modal controller of 
Arbel [2] can therefore excite vibrational modes of the 
structure. The dynamic behavior of LSS, including the actuator 
and sensor systems must therefore be carefully identified before 
deciding on the best form of the controller, its structure, and 
its gains. 

The effort here is to show the feasibility of identifying the 
modes and damping ratios of LSS by instrumental variables modal 
estimation. A simulation model of a space structure, designed by 
K. Soosaar and R. Strunce of Charles Stark Draper Laboratory 
(CSDL) (Figure 4.14) with damping provided by the low l 'Chority 
colocated controller of Ref. 2 is described by 
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where the squares of the modes of the open-loop system are 

diag Aq . [1.8, 2.77, 8.36, 8.75, 11.5, 17.7, 21.7, 
22.6, 72.9, 85.6, 106, 167] 
the control distribution matrix is 
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Figure 4.14 


CSDL Space Structure Model 
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3. 12^1 
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The system was excited by a sinusoidal frequency sweep u. The 
measurements were corrupted by white noise. Figures 4.1Sa-t show 
the transfer function measurements of the IVME algorithm, and 
Table 4.7 compares some of the actual and identified modes of the 
space structure. The dominating modes were identified with less 
than 11 error. The highly damped modes were poorly identifiable 
and resulted in larger errors up to 8% in frequency and 30% in 
damping ratio. Since the highly damped modes are naturally 
stable, they should not be considered in the reduced order model 
used for controller design. Hence, the poor behavior of the 
estimator for highly damped modes is of no practical concern. 
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APPENDIX A 

FREQUENCY DOMAIN ARMA EQUATIONS 


A linear, time- invariant dynamic system of order n is 
described by a set of linear constant coefficient differential 
equations of the form 


^ a 
i=l ®i^ 


m 


i»l 


b . u 
1 


(m-i) 


+ V 


( 1 ) 


where y is the output vector of dimension p, u is the input 
exciting the system, and the parameters a^, i=l to n and 
i=l to m are the coefficients of the characteristic 
equations for the poles and zeros of the system, respectively. 

In the frequency domain, Eq. (1) may be written as 


Y=Y ... lla . 

T 

a=[-a^, ... , 3* 
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, ... bjl’' 


u 


12 ) 


where Y and U are the Fourier transforms of y and u. 
Premultiplying each of the p equations from Eq. (2) by 


... 1]"^ , k = 1 to p and [ ( - ju)"'. . . 1 U* 

gives a set of pCp-l)/2 independent linear matrix equations 


and p equations, 

ns =S 0 ,a+SC a 11.. 
m uy^ uy^ m,n-l uu k-1 to p, 

where = Y^Y-, * U*Y^^. = U*U, 


( 4 ) 
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Equations (3) and (4) are valid for all w , and hence can be 
integrated over all w or, in the case of discrete Fourier 
transforms, averaged. The average equations are. 


Vy (f) = M (T)a + M (t)3k, k,i=l to p, 
^ki ^ki ^^k 


V » Nl , a + M 8,,, k = l to p, 
u^k “^k uu k’ 
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juJT 
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M... Ct) = I S..„ e^‘^^Q 
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‘n-l,m* 


V = in s 

u/j. u, m uXj^’ 


^‘m,n-l ’ 


M = ZS n 
uu uu m,m, 

and R is the space of all significant frequencies. 

These equations are non-singular if rank CM^jy)=m+l, and 

rank (M )= n. For the case of zero input, equation (5) 

^ki 

simplifies to 



(T)a 

^ki 


( 7 ) 


This equation defines the characteristic equation and poles of 
the system. 
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APPENDIX B 


SOLUTION OF THE INSTRUMENTAL VARIABLES 
ARMA EQUATIONS 


For a system with n poles, m zeros, p measurements and 
a scalar input u the equations defining the ARMA coefficients 
in the frequency domain are 


(t) = (T)a + (T)8y, k,i = l to p 






V * M „ ot + M B, , k = l to p, 
uyj^ uyj^ uu k ’ 


( 1 ) 

( 2 ) 


where 
of the 
delay 


V(t) and M(t) 
input and output 
( see Appendix A) . 




are defined by the auto and cross spectra 
of the system respectively. t is a time 
Solving Equation C2) for gives 
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Using .. different t, = 2nAt, results in nT*p 

=> » min ‘ 

independent equations for a . Combining these equations gives 

Aa = B (5) 
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The least 


Using 0 ( 


squares solution of Equation C5) is 
a » a'^B. 

in equation (3) gives k*l to p. 


APPENDIX C 


AVERAGING SOLUTION OF THE INSTRUMENTAL 
VARIABLES ARMA EQUATIONS 


For a system with n poles, m zeros, p measurements and 
a scalar input u the equations defining the ARMA coefficients 
and 6 in the frequency domain are: 

V (x) = M (x)ai + (x)Bj^, i = l to p, (1) 

^ki ^ki ^ 



M a + M B, , k = l to p, 
uyj^ uu k’ 


( 2 ) 


where V and M are defined by the auto and cross spectra of 
input and output of the system respectively. x is a time delay 
(see Appendix A) . 


The sum of Equations (1) is 


(t) = 

i ,k, ' ki 

Solving Equation (2) gives 
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(3) 


(4) 


(5) 


or Ca=d. The solution is a=C ^d, from which 3^ can be 
calculated using Equation (4). 
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